3 &/) / fi/y) &S 

/ Mf^L I NEAR; C ONTROL OF AIRCRAFT 
1. R. Hunt* 

Department of Mathematics 
Texas Tech University 
Lubbock, TX 79409 USA 

G. Meyer** 

NASA Ames Research Center 
MS 210-12 

Moffett Field, CA 94035 USA 
R. Su + 

Department of Electrical Engineering 
Texas Tech University 
Lubbock, TX 79409 USA 


v— 1 

C3 

TO 

• 

n 

z> m 

CO 

co 

4 

» 


r" 

cn 


03 


rn 

TO 

5 

• o 


<0 


TO 


-< 

sr 



CO 



u 

G3/08 4 

55 

UU D 

00 

NJ O 

Cfs 

<£> I-* 

1 

(1) 

NJ 

w 

00 


VJD 

Ul 

00 


ABSTRACT 

In recent work at the NASA Ames Research Center, transformations of nonlinear 
systems have been used to design automatic flight controllers for vertical and short 
take off aircraft. Under the assumption that a nonlinear system can be mapped to a 
controllable linear system, we motivate byj3arCtaT~d i f ferential equations a method 
to construct approximate transformations in cases where exact ones cannot be found. 

We also present an application of the design theory' to a rotorcraft, the UH-1H 
hel icopter. 

I. INTRODUCTION 

There are two related problems of interest in this paper. Assuming that a non- 
linear system is transformable to a controllable linear system, we first derive a 
technique for constructing an approximation of s~uch a mapping. Second, we want to 
present test results of the application of nonlinear transformations to the design 
of automatic flight controllers for aircraft. Thus the unifying thread is the 
theory of transformations of nonlinear systems to linear systems, and we shall pro- 
vide a brief review of the results in that area. 

If a nonlinear system can be transformed to a controllable linear system, then 
the transformation is presented as a solution to a system of partial differential 
equations (which can be reduced to ordinary differential equations). However, it is 
not always possible to solve these equations exactly in order to find the desired map- 
ping. By examining the partial differential equations, we derive an interesting linear 
system, called the modified tangent model , and i ndicate how an approximating transforma- 
tion (about a point in state space) can be found from this model. Results concerning 
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zfnis modified tangent model are presented in [3] arid [8]. 

r r > 

The second author first outlined his nonlinear to linear transformation • 
approach in [12], and it has been applied to several aircraft of increasing com- 
plexity. The completely automatic flight control system was first tested on a 
DHC-6 [19], and the reference trajectory for the flight test exercised a substantial 
portion of the operational envelope of the aircraft. Next, the technique was appli- 
ed to the Augmentor Wing Jet STOL Research aircraft with the successful flight tests . 
reported in [13]. Methods for providing pilbt inputs and application to an A-7 air- 
craft for carrier landing and testing in manned simulation are summarized in [15], 

[16] and [18]. Current work on the UH-1H helicopter and results are contained in 
this paper. 


II. PRELIMINARIES 


Before moving to definitions and basic theory concerning transformations of 
nonlinear systems to linear systems, we mention related references in the literature 
and the types of systems under consideration. Brunovsky [2] indicated a method for 
transforming controllable linear systems to canonical form. The second author to- • 
gether with Cicolani in [12] and [13] considered transformations of block triangular 
nonlinear systems. Krener [10] and Brockett [1] proved results concerning equiva- 
lence of systems under various type of coordinate changes. Closely related to the 
theory found in [5], [6], [7] and [17] is the work of Jakubcyzk and Respondek [9]. 


We examine nonlinear systems of the type (with controls entering linearly) 

m 

(1) x ( t ) = f(x(t)) + l u. (t)g. ( x ( t ) ) , 

i=l 1 1 

where f,g-j g are C°° vector fields on F n (or on an open set of lR n ). For our 
canonic linear system we take 


(2) y ( t ) = A Q y(t) + B n v 


0 


in Brunovsky canonical form with Kronecker indices ^ >_ 


. .>< with 
— m 


+ k 0 +. . .+<_ n. 

i 2 m 


The Lie bracket of two smooth vector .fields f and g is 
[f,g] = -| f - fig. 

where -&§■ and are Jacobian matrices. Successive Lie brackets like [f,[f,g]] = 
(ad^f,g), [f , [f , [f ,g]]] = (ad f,g), etc. can be introduced. A set of vector fields 
is i nvol uti ve if the Lie bracket of any two of them is a linear combination of the 
elements in the set. 

We are interested in the following sets, and note that the vector fields 

9 r g 2’"' ,9 m ma y be reordered. 

< 1~1 < 2-1 
C = r (ad x ,g 2 ,[f,g 2 ],. , (ad f,g 2 ) , 

<ad ™ f -V ) 


(3) 


( 4 ) 


k--2 k.-2 

*“' -P n ^ n T P n "1 f n A v 


Cj - ^>9i ) 

K • — 2 

•■.»g m »[f.g m ]....»(ad 3 f , g m ) } for j = 1 ,2,.. . ,m. 

The following local result is proved in [7]. 

Theorem 1 . The system (1) is transformable to the system (2), where the state vari- 
ables x-j .x^, . . . ,x n lie in a sufficiently small open neighborhood W of the origin in 
K n , if and only if the following three properties hold on W; 

1) the n vector fields in C are linearly independent, 
ii) the sets C. are involutive for j = l,2,...,m, and 

iii) the span of C- equals the span of C. n C for j = l,2,...,m. 

3 3 

For simplicity we assume that W = lR n , avoiding having to say for (x, ,x~ , . . . ,x ) 


K 1 


e W. We also define a new set of vector fields by adding (ad f,g n ) 
k 2 - K m 

(ad f ,g 2 ) . ,(ad f,g ) to C and reordering to obtain 


K-, -1 


{(ad ^f ,g-| ) , (ad 1 f ,g ] ) , . . . ,(ad 2 f ,g 1 , (ad 2 f,g 2 ), 


<2 - l 




(5) 


(ad f .g-j ) > (ad f ,g 2 ) , . . . ,(ad f^Mad f,g 2 ), 

K- K*3 “ 1 «o“l K^-lr- \ 

(ad f,g 3 ), (ad-- 3 f, gi ), (ad f,g 2 ), (ad 3 f ’V’ 

- - . .g-, »g 2 »- • • ’9 m }- 


m 


Before 


e forming this set we should check k, > k 0 >...>k . If k, = k„, then 
■, I — 2— — m i 2 

“ I <2 Kj 

(ad f.g^) appears after (ad f,g 2 ) = ( ad f,g 2 )> etc. 

For the remainder of this paper we assume that system (1) satisfies the hy- 

potheses of Theorem 1. Next we consider the problem of constructing transformations. 

III. MODIFIED TANGENT MODEL 

To build a transformation taking system (1) to system (2) we must solve the 

partial differential equations (here a-, = tc-, , a 0 = k,+k 0 , a 0 = k,+k 0 +k_ , . . . ,a = 

r I I 2 I 2 J J . 2 3 m 

m 


k 1 +k 2 + " ‘ - +Km ” ’ 


<dT^, (ad"*f,g^)> = 0,j = 0,1,... ,*^-2 and i = 1,2,. ..,m, 
<dT • ; (ad^f ; 9^ ) ^ 0 , g 0,1, .»»,£*■— 2 and i — 1,2,.*., in, 

O » tX jl ^ 


<dT 


a +1 ,(ad J f ,g i )> = 0, j = 0,1, . . . ,^-2 and i = 1,2, ... ,m, 
m-l 


(6) <dT^ , f > — T^__^ , k ~ 1,2,.. . ,a ; k ft fOj / • • • , , 

m -1 

<d T , f > + j V dT l' (aa £ '9 i > > = T n+1- 
1 1=1 


m 


< 2 -i 


<dT0 2 ,f> ± (ad f, 9i )> = T n+2( 


m k -1 

<dT n' f> i . I .V a, o 1+ l' (adm f '3i )>=T n +m 

1 = 1 m- l 


with the, matrix 


(7) 


Kl-i 


K r x 


<dT^, (ad f , g^) > . . . <dT^ , (ad f,g m )> 
<aT o 1 + i ' <ad2 (ad 2 f,g m )> 


<dT a 1+ l' (adI ” f >9 1 )>...<dT + i-< adnl f.S„>> 

m- 1 m-l 


being nonsingular. Here + is for ^ odd and - is for k. even. The important 
equations in (6) to consider are the first. m, the others following by easy Lie 
(differentiation. Of course, it is not always possible to solve the ordinary 
differential equations in closed form, and that is the reason for the method to be 
introduced shortly. We want to emphasize a related set of partial differential 
equations that a transformation T must solve. These can be deduced from (5) and a 
Leibniz formula (see [7]). 


<dT 1 , (ad^f ,g i ) > = 0, 

j = 

0,1,. ..,<^-2 and i = 

1 ^ 2 ^ * • • f xn f 

<dT 2 , (ad-^f ,g i ) > = 0, 

j = 

0,1, . . . ,ic^-3 and i = 

1,2,. • . f xn f 

<dT a 1 -l'5 i > = 0, i = 

1,2, 

* • • ,rn , 


<dT o 1+ l'< adjf 'Si>> = 

0, j 

0,1,... , k 2 ~2 and x 

1,2,# •«, xn , 

(8)<dT a +2 , = 

l • 

• 

• 

0, j 

= 0,1,..., < 2 ~3 and i 

1^2^##«j xn , 

<dT o 2 -i'9 1 > ; i = 

1,2, 

. . . , m 


<dT o , +1 .(ad :l f,g i )> 

= o, 

j 0 , 1 , . . . , ic “2 and 

i - 1 , 2 , * . • ,ra 



dT 0 Jn . 1 +2' (aa; ' f '9i) > = 0, j = 0,1 < m -3 and i = 1,2 


<dT 0 -l'^* ~ i - 1,2, ...,m, 
in 

and we again want the matrix (7) to be nonsingular. 

With these equations in mind we turn to our development of the system called 
the modified tangent model. A comparison of this model and the tangent model [11] 
will be made later in this section. 

Suppose we take a pointx Q in state space and linearize a transformation T about 
x Q and denote the linear part by T x . Then solves equations (8) where each 
(ad^f,g.j) is replaced by (ad^f ^ ) (x q ) . The motivation behind the modified tangent 
model is to find a linear system 

(9) x = f (x ) - Ax + Ax + Bu 
o o — 

so that the Lie brackets of Ax and B= ( by ; b^ » .' '.b m ) agree with the corresponding 

brackets of f and g. ,g 0 g at x . 

3 1 3 m o 

Definition 2 . The system (9) is called the modified tangent model at x q for system *■ 
(1) if A and B satisfy (take + for k even and - for k odd) 


A b ] = +(ad k f ,g-j )(x Q ) , k = 0,1,. 

A k l> 2 = +(ad k f,g 2 )(x o ) , k = 0,1 ,..., k 2 


A k b m = ±(ad k f,g m )(x 0 ) , k = 0,1,..., K m 

A transformation taking Ax + Bu into Brunovsky .canonical form (2) satisfies 
equations (8) where each (ad J f,g i - ) is evaluated at x . 

There is an interesting geometrical interpretation of the modified tangent 
model. Letting C(x q ) and Cj(x Q ) denote the evaluations of the sets in (3) and (4) 
at x Q , we find that for every jJ<j<m, C ^ ( x o ) n C(x q ) spans the tangent space at x q 
of the integral manifold of C; n C guaranteed by the Frobenius Theorem. Since 

ki -1 V 1 

C(x ) = (b 1 ,-Ab 1 ,. . . ,+A 1 b ] ,b 2 ,-Ab 1 . . ,+A b 2> .... 


K -1 

b ,-Ab , ... ,+A m b m } 
mm— m 


k --2 k --2 k ^-2 

C i X 0 ) = b i >“‘ Ab 1 ,. . - ,+A 3 b 1 ,b 2 ,-Ab 2 ,...,+A J b 2 ,...,b m ,-Ab m ,. . .+A b m ) 

j— 1,2,... ,m. . . 

the same is true for the correspondi ng brackets generated by equations (8). 

We remark that it is extremely easy to calculate A and B from (10) despite the 
nonlinearities of the equations involved. We form the (n+m)(n+m) matrix whose first 
column is the first vector in (5) plus m. zeros, second column is the second vector 
in (5) plus in zeros, etc. 



*--ijse orthogonal coordinate changes on ]R n to take this matrix to "generalized 

lower Hessenberg form" as in [4]. The linear independence of b m ,b m -j . . . : ,b^ .Ab^, 

^m-l »• •• » Ab-j • »A ^ b^ (i , e . the vectors in C(x )) allow us to trivially solve 

in order for b ,b , ,...b., the last column of A, the (n-l)th column of A)..., the 
m m- 1 I 

first column of A. 

The purpose of the theoretical part of this paper is to construct an approximate 
transformation taking equations (1) to equations (2). This approximate mapping is 
found by giving appropriate solutions to equations (8), where each Lie bracket is 
evaluated at x . As remarked earlier, the transformation that moves Ax+Bu in the 
modified tangent model to Brunovsky form is designed to solve those equations. A 
simple algorithm for taking a controllable linear system to canonical form is pro- 
vided in [4]. It involves taking the modified tangent model to a "block triangular 
system" (see [4]), and the orthogonal coordinate changes mentioned earlier accomp- 
lish this. Next,' "Lie differentiation" completes the linear transformation for 
Ax+Bu. For the approximate transformation, this Lie differentiation is done with 
respect to f(x Q )-A(x )+Ax. 

The tangent model approach was introduced in [11] and applied on a trajectory 
autopilot for a helicopter. For the tangent model we write equations in the same 
form as (9) 

(11) x = f (x )-Ax +Ax+Bu, 

0 0 

except that A is the Jacobian matrix of f(x) at x Q , and B=(b 1 , . . . b^) is an nxm 
matrix with 9 i(x q ) as its first column, g^x^ as its second column, ... »9 m ( x 0 ) as its 
last column. The first step in constructing an approximate transformation by use of 

(11) is to take Ax+Bu to linear canonical form. However, it is not in general true 
that -[f,g^](x o )=Ab-j, [f ,g^](x 0 )=Ab 2 , etc. (an example is provided in [3]). Hence, a 
transformation derived in this manner may not have the desirable geometric inter- 
pretation possessed by the modified tangent model. However, the tangent model seems 
to provide an adequate transformation in many cases, as illustrated by its use on 
the simulator for a tail-sitter aircraft. ' 

IV. APPLICATION'™ A HELICOPTER 

The helicopter will be represented by a rigid body moving in three-dimensional 
space in response to gravity, aerodynamics , and propulsion. The state, 

(12) x = (r,v,C,w) T e X c JR 3 x 1R 3 x S0(3) x 1R 5 

where r and v are the inertial coordinates of body center-of-mass position and veloc- 
ity, respectively, and C is the direction cosine matrix of the body-fixed axes rela- 
tive to the runway-fixed axes (taken to be inertial). The attitude C moves on S0(3). 
The body coordinates of angular velocity are represented by w. 


The controls are 


M 

where u is the three-axis moment control, that is, roll cyclic and p-itch cyclic, 
which tilt the main-rotor thrust, and the tail-rotor collective, which controls the 

p 

yaw moment; and u is the main-rotor collective, which controls the main-rotor thrust. 

The effectively 12-dimensional state equation consists of the translational and 
rotational kinematic and dynamic equations: 
r = v 

v = f F (x,u) 

(14) . 

C = S(wjC 

O) = f (x ,u) 

F M 

where f and f are the total force and moment generation processes, and (x,u) are 
defined by Eqs. (12) and (13). 

M 

In general, the moment generation process f is invertible with respect to the 
M 

pair (uj,u ), and for the restricted class of maneuvers being considered in this 
experiment (i.e., no 90 rolls), f is invertible with respect to the pair(v 3 ,u ). 

So the four controls u may be related to four accelerations (w,^) by a nonsingular 
transformation, say, 

(15) u = h M (x',(<L,v 3 )) 

If (to, Vo) a re chosen to be the new independent control variables (a,aj to replace 

•j m p J 

the natural controls (u ,u ), then the state equation (14) becomes the following: 

r = v 

I ; i \ , 

= f°(r,v,C) + f [r,v,C,(o,(a,a J] 

\ V 2 

(15) v 3 = a 3 

C = S(to)C 
<0 = a 

The function f^ in (15) represents (nonlinear) zeros: 

(17) f 1 = (r,v,C)u> + £ 2 (r,v,C)a + e 3 (r,v,C)a 3 + f 2 

o 

where f is second or higher order in (w,a,ao)- In the case of the actual heli- 
2 J 

copter, f is very small and completely negligible. Consequently, the first condition 
for linearizability, which requires that controls (a,a 3 ) enter linearly into the state 
equation is in fact satisfied. With the new controls, the state equation becomes, 



r = v 


( 18 ) 


f v, \ 

I • J= f°(r,v,C) + e-j (r,v,C)(u + e 2 (r,v,C)ct + e 3 (r,v,C)a 3 

v 3 = a 3 
C = S(u)C 

to = a 


The primary means for controlling the horizontal motion of the helicopter is through 
f°. The e. teams are parasitic and negligible for regulator bandwidths below 0.5 
rad/sec. The dominant team f° is invertible with respect to the pair ((v-j v 2 ,E 3 ),C) 
in which E 3 , a rotation (angle \p) about the veritical, defines the heading of the 
helicopter. So, the hel i copter, attitude C can be related to horizontal acceleration 
and heading (a^,a 3 ,E 3 ) by a nonsingular transformation, say, 

(19) C = h F (r,v,( ai ,a 2 ,E 3 )). 

From (18) with £^.=0 and (19) it follows that angular velocity 

F T * 

(20) u = h u (r,v,a,E 3 )(a 1 ,a 2? ji) 

where h^ is nonsingular and a=(a^ ,a 2 ,a 3 ) T , T now denoting transpose. 

In the actual implementation of the control scheme, the inverse transformations 
h M and h ^ do not appear explicitly. Instead the Newton-Raphson algorithm is used 
to compute u from f n and C from r . The Jacobian matrix h w in (19) is available as 
a by-product of the algorithm. 


The canonic model has the Kronecker index set (4, 4, 2, 2}. 

5 

y and control y are chosen to be the following 


(21) 


y 1 = (r r r 2 ) T 

y 2 = ( V v 2 ) t 

3 T 

y — (^-j 

4 r • -vT 
y ^i*^2* v 3 * ^ 

y = (a 1 ,a 2 ,a 3 ,^p) 


The canonic coordinates 


5 

The map from canonic controls y back to natural controls u is given in two steps. 
First, commanded angular acceleration is computed by taking the time derivative of 

• p 

(20) along the model trajectory, y(t), and neglecting h , 


(22) a = h^(r,v,a,E 3 )(a 1 ,a 2 ,^) T 


Then the controls u are computed from (15), 
(23) u = h M (x ,(a,a 3 )) 

where x is the natural state defined by (12). 



Extensive testing indicates that, at least in the present application of the 
particular helicopter, the two approximations f^=0 and li^rO, which greatly reduce 
the computing load on the flight computer, cause insignificant errors.' The design is 
robust with respect to modeling errors such as weight and location of the center of 
mass. However, in high speed flight there is too much sensitivity to errors in wind 
estimates. The problem is due to the independent control of heading, which, in high 
speed, should point the helicopter in the direction of the relative air velocity 
vector. The uncontrolled helicopter has a natural tendency to point into the total 
wind. This and other benign built-in characteristics are removed by the transformation 
to Brunovsky form. A control system structure that allows the retention of the 
benign characteristics of the plant is shown in Figure 1. 



Figure 1. Structure of the Generalized 
Model Follower. 


The transformations T -|3’ T 23 ,W 13 ,W 23 are nons ' in 9 u1ar - The usual configuration 
(see [13]) is then a special case in which the model dynamics u^-^x™ and regulator 
dynamics u™— >x™ are both Brunovsky (i.e. T^ and are identities). With the 
additional freedom provided by the new structure, it was possible to retain the 
benign weather cock stability of the helicopter and thereby reduce the system 
sensitivity to wind estimator errors to tolerable level. 
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